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THE VARIABLE WARM ABSORBER IN CIRCINUS X-1 
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ABSTRACT 

We observed Circinus X-1 twice during a newly reached low- flux phase near zero orbital phase using 
the High-Energy Transmission Grating Spectrometer (HETGS) onboard Chandra . In both observations 
the source did not show the P Cygni lines we observed during the high-flux phases of the source in 
2000 and 2001. During pre-zero phase the source did not exhibit significant variability and exhibited 
an emission-line spectrum rich in H- and He-like lines from high Z elements such as Si, S, Ar, and 
Ca. The light curve in the post-dip observation showed quiescent and flaring episodes. Only in these 
flaring episodes is the source luminosity significantly higher than observed during the pre-zero phase. 
We analyzed all high resolution X-ray spectra by fitting photoionization and absorption models from 
the most recent version of the XSTAR code. The pre-zero phase spectrum could be fully modeled with 
a very hot photoionized plasma with an ionization parameter of log ^ — 3.0, down from log ^ = 4.0 in 
the high-flux state. The ionization balances we measure from the spectra during the post-zero phase 
episodes are significantly different. Both episodes feature absorbers with variable high columns, ionization 
parameter, and luminosity. While cold absorption remains at levels quite similar to the one observed 
in previous years, the new observations show unprecedented levels of variable warm absorption. The 
line emissivities also indicate that the observed low source luminosity is inconsistent with a static hot 
accretion disk corona (ADC), an effect that seems common to other near edge-on ADC sources as well. 
We conclude that unless there exists some means of coronal heating other than X-rays, the true source 
luminosity is likely much higher and we observe obscuration in analogy to the extragalactic Seyfert 
II sources. We discuss possible consequences and relate cold, luke-warm, warm, and hot absorbers to 
dynamic accretion scenarios. 

Subject headings: stars: individual (Cir X-1) — stars: neutron — X-rays: stars — binaries: close — 
accretion: accretion disks — techniques: spectroscopic 



1. INTRODUCTION 

Not many X-ray binaries show as large a variety of X- 
ray featur es and properties a s Circinus X-1 . Since its 
discovery (jMargon et al.|[T971f ). this X-ray binary has ap- 
peared at various X-ray brightness leve ls ranging from 
only a few to several hundreds of mCrab (jParkinson et al.l 
l2003f ). Its long-term brightness featured flaring activ- 
ity at increasing levels during the 1970s and 1980s and 
a steady increase in overall brightness during the 1990s 
until at some point it reached a flux of about 1.5 Crab. 
The discovery o f broad P Cygni Une profiles during these 
brig ht periods (|Brandt fc Schuld l2000l : ISchulz fc Brandt) 
l2002i) demonstrated the presence of a high- velocity outflow 
likely in the form of an accretion disk wind . The simul- 
tanuous pres ence of possibly relativistic jets (jFender et al.l 
Il998l l2004af) underlines the highly violent nature of the 
source at the time. 

Since then Cir X-1 has calmed considerably in its X- 
ray emission and although the X-ray light curve still shows 
its characterist ic pattern wit h an orbital period of 16.6 
days (fKaluzien ski et al.lll976f ). its outbursting power near 
zero phase also has diminished to only a fraction of the 
levels it showed a few years ago. In recent months the 
source was only barely or not detected with the RXTE 



ASM (jJonker et al.l [2007af ). In most repects the nature 
of Cir X-1 is still poorly understood. Despite advances 
in recent years, there remains uncertainty about even the 
most basic properties of this system. Recent observations 
of quasi-periodic oscillations (QPOs) added more incen- 
tives to a long list of indirect evidence that the compact 



object in Cir X-1 is in fact a neutron star ( Shirev et al 



Tauris et al 



1999; Tc nnant et a l. 1986: j Qu et"all 120011 : 
199i^ iBoutloukos et al.n2006[ ). Some new insights on the 



nature of its compan ion came from a recent analysis by 
iJonker et al.l (l2007br) who argue for an earlier A to late 
B-type star. iMonetil (|19920 had detected a heavily red- 
dened optical counter part with strong, a s ymmetric and 
variable Ha emission (WhelanetaL l ll977t (M ignani et al.l 
[r997; Johnston ct al. 2001). 

Th e system's orienta tion has become a much debated is- 
sue (iSchuIz et al.ll2006f) . Perhaps one of the most exciting 
but also controversial features of Cir X-1 is the existence of 
a possibly ultrarelativistic radio jet. Under the assumption 
that the jet was launched during the most recent perias- 
tron passage, iFender et all (|2004a[) deduce an upper limit 
to the inclination of about 11°. Most recent projections 
of or bital parameters pos tulate an ang le of no le s s than 
14° (! Jonker et all l2007bl ). Recentlv Heinz et al.l (|2007D 
found evidence for a parsec scale X-ray jet from Cir X- 
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1 which also appears consistent with some lower inclina- 
tion for an unbended jet. X-ray spectral studies do not di- 
rectly provide evidence about inclination since no eclipses 
are observed. However, extensive studies with ASCA and 



RXTE (iBrandt et al.|[T9 96: Shir ev et al lll999D. and later 



with Chandra (Brandt fc Schulz 



y et au\iy\)\)[]. and later 
120001: ISchulz fc Brandt! 



|2002| ) produced arguments for a highly inclined viewing 
angle including an equatorial wind-like outflow that now 
seems to be a quite typical property of these types of X-ray 
binaries (|Blundell et al.l[200lHMriler et al.l[2006allbl ). 

In this respect Circinu s X-1 bel o ngs t o a now long list 
of Galactic microquasars iMirabell (|200lh named in anal- 
ogy to their extragalactic cousins, the Seyfert galaxies. 
One of the trademarks of active galactic nuclei is, be- 
sides an X-ray bright nucleus and prominent relativistic 
jets, the existence of warm absorbers. There are vari- 
ous types of X-ray absorbers in X-ray sources. Widely 
known in Galactic X-ray binaries is photoelectric absorp- 



lar medium (iMorrison fc McCammon 


19831 Wilms et al.l 


120001 iPaerels et al.l 20011: lJuett et al.l 


2004', '20061 as well 



120031: ISchulz et all l2002al : lJuett et all I2006D . The most 
prominent indication for the existence of large amount of 
such cool circumstellar material is the existence of signif- 
icant line fluorescence as well orbitally correlated column 
densities. The existence of Seyfert-like warm absorber out- 
flows in microquasars other than Cir X-1 has recently also 
bee n observed with Ch andra in GX 339-4 and GROJ 1655- 
40 (jMiller et al.ll2004 [2b06b). The term warm absorber 
was introduced during the mid-1980s bv iHalpernl (|1984( ) 
as an illuminated plasma of moderate electron tempera- 
ture (< 10^ K) which shows similar levels of ionization as 
a collisionally ionized hot gas. Quite impressive examples 
of war m absorbers hav e been shown for AGNs like MCG6 - 
30-15 (iLee et al.ll200l) and NGC 3783 (|Kaspi et al.ll2002fl . 
A recent over view of warm absor ber properties in AGN 
was given by iBlustin et all (|2005D . 

In this third paper, we present results from two new ob- 
servations of the Galactic microquasar Circinus X-1 per- 
formed during a phase where the source was at an over- 
all extremely low flux level. A s in papers I and II 
(|Brandt fc Schuld l2000l (BSI) and ISchulz fc Brandtl[200^ 
(SBII)) the observations were performed with the High- 
Energy Transmission Grating Spectrometer (HETGS) on- 
board Chandra at "zero phase." Zero phase in Cir X-1 is 
thought to be associated with the periastron passage of the 
neutron star. The observations include a shorter exposure 
during the pre-zero phase dip and a longer exposure at zero 
phase where the source usually exhibits strong variability. 
The second of these observations was also t he basis for 
the discovery of the parsec-scale X-ray jet bv lHeinz et al.l 
12007). Througho ut this paper, we adopt a dista.nce to 
Ch X-1 of 6 kpc (jStewart et al.lll993t ICase et alJ [19981) 
and an interstellar column den sity of about 2 x 10 
Ce.g.. lPredehl fc Schmittl [19951) . 
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Fig. 1. — The light curves of obs. Ill (top) and IV (bottom). The 
data were binned by adding all HETG 1st order photons from 1.6 
to 10 A. The y-scale of the two curves was kept the same for better 
comparison. The inset in the top image shows a blow up of III. 



2. CHANDRA OBSERVATIONS 

Cir X-1 was o bserved with the HETGS (see 
ICanizares et al.l [20051 for a detailed description) on 2005 
January 3 (starting at 01:14:52 UT) for 25 ks (observa- 
tion III, OBSID 6148) and on 2005 June 2 (starting at 
23:44:51 UT) for another 50 ks (observation IV, OBSID 
5478). The denominations 'IIP and 'IV' stem from the 
fact that we observed Cir X-1 a few years earlier in a very 
similar manner and instrument configuration which were 
labeled as observations I and II (SBII). Both observations 
were performed near zero phase of the binary orbit. Ob- 
servation HI was performed during an often occurring pre- 
zero phase dip at phase 0.96 - 0.99, while obs. IV occurred 
a few phase bins la.ter at phases 1.06 - 1.10 (adopting the 
ephemeris of [Glasi[l99^. Both observations were config- 
ured the same way. The detector array (ACIS-S) consists 
of 6 CCDs which were clocked choosing a 512 pixel row 
subarray placed near the readout in order to reduce the 
frame time to about 1.7 s. There were two distinct differ- 
ences in the setup of the two observations which stemmed 
from two erroneous settings during obs. HI. While the first 
produced a different off-axis angle of the telescope with the 
only consequence that chip gaps in the effective area fall 
into different wavelength bands, the second was a bit more 
costly. Due to an incorrectly placed exclusive count win- 
dow filter on the central CCD (and we detect no source 
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Obsid 


label 


Start Date 


Start Time 


Exposure 


HEG -1st rate 


Phase 






[UT] 


[UT] 


[ks] 


cts s^-'- 


range 


6148 


HI 


Jan 03 2005 


01:14:52 


24.48 


0.228 


0.96 - 0.99 


5478 


IV 


Jun 02 2005 


23:44:51 


51.62 


0.524 


1.06 - 1.10 



photons above 10 A), the spectra in obs. HI possess only 
about 50% of the anticipated exposure. 

All observations were reprocessed using CIA03.4 with 
the most recent CALDB products. In order to fine tune 
the wavelength scale we redetermined the zero-order po- 
sition in both observations. For obs. IV we have a strong 
0th order point source with read out trace and positive and 
negative order dispersion. This allows us to reach a posi- 
tional accuracy of about 1/4 of a detector pixel ensuring 
a wavelength scale accuracy of about 1/8 of a resolution 
element, i.e. 0.0025 A for MEG and 0.0015 A for HEG 
spectra. In obs. HI the 0th order was filtered for every 
10th photon and weak, and due to the misplaced exclusion 
window we do not have a read out trace and significant 
negative order spectra. This severely limited our ability 
to determine a precise 0th order position. Therefore for 
obs. Ill we only reach an accuracy of 0.010 A for MEG and 
0.005 A for HEG in our wavelength scale. For transmission 
gratings the dispersion scale is linear in wavelength and we 
perform the spectral analysis always in wavelength space 
to avoid non-linear binning. We used standard wavelength 
redistribution matrix files (RMF) but generated ancillary 
response files (ARFs) using the provided aspect solutions, 
bad pixel maps, and CCD window filters, [j For all the 
observations we generated spectra and analysis products 
for the 1st orders only. 

Figure [T] shows the light curves of the two observations. 
Note that since obs. HI (top) occurred at pre-zero orbital 
phase, the count rate appears exceptionally low. However, 
even compared with observation 11 from the year 2000 at 
similar phase, the rate has dropped by over two orders 
of magnitude. Observation IV (bottom) shows not only 
the expected higher rate, but also quite significant vari- 
ability. Besides episodes of relative quiescence the source 
now engages in frequent flares. Though this is not unex- 
pected during this orbital phase, the observed activity has 
a somewhat different appearance compared to the activ- 
ity in previous years where the fiux rose rather quickly 
to a high level accomp anied by many dipping episodes 
(e.g.. IShirev et al.lll999f ). Compared to the fiux levels in 
previous years, this rate of the source even in this phase is 
now quite moderate and also reduced by almost two orders 
of magnitude. 

2.1. Continuum Spectra and Fluxes 

In order to determine the source flux during the two 
exposures we fit the source spectra with generic contin- 
uum models. Unlike in obs. 1 and II, where we had severe 
spectral degradation due to photon pileup, we do not have 
to restrict ourselves to specific spectral ranges or rely on 

^see |http: //asc .harvard.edu/ciao/threads/ | 



higher orders in our new observations. While obs. HI is 
practically pileup free, possible pileup at the maximum 
photon fiuxes in obs. IV is below 2%, which we treat as 
a contribution to the systematic uncertainty. For the flux 
analysis we also exclude detected line features (see below) 
which although not significantly contributing to the over- 
all source fiux will affect the broad band continuum fit. 
The main motivation for this portion of the analysis is to 
obtain X-ray fluxes which serve as a check on the outcome 
of the fits using the detailed photoionization model below. 
We also investigate potential similarities of the continuum 
shape with previous observations. 

Figure [2] indicates that there is no significant fiux long- 
ward of 10 A. This is different from the previous obser- 
vations, where we have flux contributions up to 14 A and 
use a model sp ectrum based on previous fits to ASCA data 
dBrandt et al.l 1996) consisting of two blackbody compo- 
nents with both intrinsic partial-covering absorption in 
addition to interstellar absorption. The raw spectra in 
obs. Ill and IV indicate a high level of complexity involv- 
ing a combination of emission and absorption. We could 
not produce good fits for both observations using the dou- 
ble partially covered blackbody model, but found similar 
options. Table 2 summarizes the fit results. Statistically 
only one simple power law plus absorption was required for 
observation HI. This was not the case in obs. IV, which in 
order to include the softer end of the spectrum (i.e., above 
5.2 A ) needed an additional component. iGallowav et al.l 
(|2007f) successfully fitted a single blackbody with partial 
coverage absorption to high fiux spectra at higher orbital 
phases. The fits produced reduced below 1.5 only under 
certain conditions. In obs. HI the fit only allowed fdor in- 
terstellar absorption leaving significant residuals near the 
Si and Fe K edges. In observation IV the model fit only if 
the wavelength range was restricted to 2 - 8 A. 

The next step involved a two-component model consist- 
ing of two power laws, also with additional absorption. 
One of the components required significant additional ab- 
sorption with a partial covering fraction above the 90% 
level. Substitutions of one of the components with black- 
body functions left the reduced well above 1. Still, the 
two-component configuration is relatively similar to what 
we found in observation I and II (SBII) and we there- 
fore applied it to obs. HI as well. We found a good fit 
with rather similar parameters. Here the substitution of 
the absorbed power law with a blackbody worked as well. 
The conclusion we draw from this is that the new spec- 
tral shape is likely not much different from the one in the 
high-flux states a few years ago. The models also sug- 
gest that the transition between pre- and post-zero phases 
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Fig. 2. — The raw spectra of obs. Ill (top) and IV (bottom) from all HETG 1st order photons. In both cases the spectral bandpass is 
limited to wavelengths below 10 A indicating high column densities. 
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is predominantly a change in the hard absorbed spectral 
component. Most important, however, is the fact that 
the continuum is consistent with a power law shape in 
both, the hard and softer parts of the spectra, which then 
allows a relatively straightforward treatment of photo- 
ionization and -absorption. The absorbed source flux (2 - 
10 keV) in obs. Ill was determined to (1.23±0.02) x 10-^° 
erg cm~2 and in obs. IV was (5.72±0.03) x 10-^° 
erg cm~^ s^^ . The interstellar column density is assumed 
to be 1.94x 10^^ cm~^ as provided by many previous X-ray 
analyses and which is close to the pre dicted column from 
H I m easurements in this line of sight (jPredehl fc SchmTtl 
Il995l) . Both spectra in all models show that substantial 
additional absorption is required. 

3. DETAILED SPECTRAL ANALYSIS 

Our models for the spectral features (emission and 
absorption lines, and bound- free absorption) consist of 
both phenomenological fits, using Gaussian lines, and fits 
to photoionization models constructed using the XSTAR 
(v. 2.11n2) code0 XSTAR calculates both the ioniza- 
tion/excitation of a gas illuminated by a strong source of 
X-rays and also the spectral features imprinted by the gas, 
both in emission and absorption. This is available as an 
'analytic' model in XSPECQ, which can also be imported 
into ISIS|!| The models assume a time-steady balance be- 
tween processes affecting ionization, recombination, exci- 
tation, de-excitation, heating and cooling of the gas. The 
key parameters determining the emission and absorption 
properties are the shape of the incident spectrum and the 
ionization parameter, defined as ^ = Airfx/ue, where X- 
ray flux is defined as the ener gy flux incident on the 
gas integrated from 1 - 1000 Ry (jTarter et al.lll969f ) and 
rig is the electron number density in the gas. Additional 
parameters include the elemental abundances, bulk veloc- 
ity or redshift, and turbulent velocity, in addition to the 
column density or the emission measure for absorption or 
emission models, respectively. Each model component is 
idealized as a uniform medium with a given set of input 
parameters. Emission and absorption components can be 
combined, either singly or multiply, using standard 'model' 
commands in XSPEC or ISIS. 

For all spectra we performed simultanous fits to the in- 
dividual 1st order spectra of the MEG and HEG. Only 
for plotting purposes are the spectra and model fits added 
together. During the analysis we binned the spectra into 
0.005 A bins which matches the available MEG bins and 
undersamples the HEG by a factor of 2. 

3.1. Emission Lines 

Table 3 lists all detected lines with positions and fluxes 
obtained from Gaussian fits. The values were obtained by 
fixing the continuum level to the fits from Sect ion [TT] The 
line lists for the two observations differ substantially due 
to the fact that in obs. IV we do not find any line emis- 
sion below 3.9 A except for traces of an Fe K fluorescence 
line. In obs. HI we detect all major H- and He-like Lyman 
a lines plus a few weak H-like Lyman /3 lines from some 
elements. The fits to both spectra indicate no systcm- 

® http: / /hcasarc. nasa. gov /docs/ software / xstar / 
'^http: / /heasarc. nasa. gov /docs/ software / xspec / 
8http;//space.mit.edu/CXC/ISIS/ 



atic or significant line shifts. Observation HI appears to 
be slightly inconsistent with this assessment. The labora- 
tory wavelengths tabulated in column 2 are taken from ta- 
bles in the Atomdb web-guide. The difference in measured 
wavelength from the laboratory values for observation HI 
distribute fairly well around a centroid shift of 0.0045 A, 
which is well within the statistical and systematic uncer- 
tainties expected for this data set. In obs. IV there are 
fewer lines but they distribute well around the expected 
rest wavelengths. 

The line shapes are resolved but require little broad- 
ening beyond the resolution of the HEG. Most Doppler 
broadening line velocities are less than 800 km s~^, except 
for Fe and some Mg lines. The 90% confidence limits in- 
dicate that the lines themselves are just barely resolved. 
There is no obvious correlation with wavelength scale ex- 
cept for the Fe region where the resolution is already over 
1500 km s~^. The values for Mg stand out a bit. Specifi- 
cally the fact that we could not resolve the Mg XI triplet 
indicates some higher velocity broadening. Since in both 
observations the lines are just barely resolved it is also not 
obvious if the trend to lower velocities in obs. IV is statisti- 
cally relevant. Although the lines tend to be brighter, the 
fits suffer from lower contrast due to a higher continuum. 

Another peculiarity between the line emission of obs. HI 
and IV is the fact that we do not detect any emission lines 
with wavelengths shorter than Ar XVII in observation IV. 
Using the 1 cr statistical uncertainty of the underlysing 
continuum we limit the line fluxes for the undetected line 
top below 1.2x10"'' ph s"i cm . While we do not detect 
Ar XVIII and Ca XIX in either absorption or emission, 
we detect Ca XX, Fe XXV, and Fe XXVI in absorption. 
Furthermore there is an unresolved line absorption com- 
plex longward of Fe XXV, likely due to lower transitions 
of Fe XX to XXII indicating the presence of both highly 
ionized and moderately ionized absorbers in the line of 
sight. We refer to these respectively as 'hot' and 'warm' 
components, although these designations do not reflect the 
fact that the dominant source of ionization is likely the 
continuum photons from the compact object rather than 
thermal electrons. The line positions of the hot absorp- 
tion lines do not deviate significantly from the laboratory 
position. This is different to obs. I and II where they ap- 
peared blue-shifted, thereby producing P Cygni lines. A 
more detailed analysis of the absorbers is given below. 

The Fe K region shows some peculiarities as well. In ob- 
servation HI the Fe XXV is clearly the dominant feature. 
It does show some excess fiux longward of the forbidden 
line position, however with low statistical significance. Un- 
der the assumption that we are observing a photoionized 
plasma, then these could be trace contrib utions from in- 
ner shell transitions of Fe XIX - XXII ions (jKallman et al.l 
2004). If so, this would suggest the presence of a distri- 
bution of ionization parameters, extending to values ^10 
times less than the dominant component. 

The main challenge for the photoionization model comes 
from the He-like lines, here Fe XXV, Ca XIX, Ar XVII, 
S XV and Si XIII and Mg XI, observed in obs HI. In in- 
terpreting these lines we noted their dependence on gas 
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TABLE 2: THE CONTINUUM FIT RESULTS 



Obsid Model 



[10^^ 



cm 



cm 



[cts 



Al 

1- 



r2 or kTbb 



[cts 



A2 

s-1 A- 



6148 P 

PCBP 
PCPP 

5478 PCPP 



o O«+0.30 
0.00_Q 29 

^•^^-0.08 



+0.45 
0.64 
+0.08 
-0.05 



4.46 
1.84 

12.67t°-;? 



0.95iro^ 
0.92iro^ 
0.95tHl 



2.11 



+0.16 



-0.15 
1+0.03 
-0.16 
hO.03 
-0.01 

1 67+°-°° 
■'■•°'-o.oi 



1.40] 
1.87j 



2.3tEJ;^ X 10" 



6.8 
1.6 



-0.2 
-1.8 
-0.6 
-0.1 

n 26+°°° 



X lO"'* 

X 10^2 



7.08t?:°^ 
5.99t°;i 



6.lt^:« X 10-3 
5.3t^ ° X 10-3 
8.4;?:° X 10-2 



1.14 
0.90 
0.94 
0.84 



P = Power law; PCBP = partially covered blackbody + power law ; PCPP = partially covered power law + power law 



TABLE 3: X-RAY LINE PROPERTIES 



Observation III 



Observation IV 



ion 


Ao 




Am 


Fluxiii 


Vlll 


Aiv 


Fluxiv 






A 




A 


10-^ ph s-i cm-2 


km 


A 


lO-'' ph s-i cm-2 


km 


Fe XXVI 


1.780 


1, 


.783±0.002 


0.97±0.20 


1680±1490 


1.776±0.001 


-2.45±0.08 


740±880 


Fe XXVr 


1.850 


1, 


.853±0.001 


1.58±0.17 


570±1130 


1.848±0.002 


-1.78±0.34 


1090±850 


Fe XXVi 


1.859 


1, 


.860±0.002 


0.80±0.15 


1520±1200 








Fe XXVf 


1.868 


1. 


.870±0.001 


1.07±0.15 


570±1120 








Ca XX La 


3.021 


3 


.023±0.002 


0.26±0.06 


210±450 


3.019±0.001 


-0.83±0.19 


430±520 


Ca XlXr 


3.177 


3 


.163±0.012 


<0.07 


<650 








Ca XlXi 


3.189 


3 


.184±0.001 


0.51±0.15 


660±950 








Ca XlXf 


3.211 


3 


.205±0.003 


0.25±0.10 


590±650 








Ar XVIII La 


3.734 


3, 


.738±0.005 


0.40±0.13 


590±990 








S XVI L/3 


3.780 


3 


.787±0.006 


0.24±0.12 


550±470 








Ar XVIIr 


3.949 


3 


.939±0.004 


<0.06 


<610 


3.932±0.005 


<0.01 


<660 


Ar XVIIi 


3.966 


3, 


.956±0.004 


0.24±0.13 


360±530 


3.958±0.007 


0.19±0.18 


130±110 


Ar XVIIf 


3.994 


3 


.995±0.003 


0.45±0.15 


600±520 


3.998±0.005 


0.34±0.19 


600±490 


S XV L/3 


4.310 


4, 


.310±0.006 


0.28±0.15 


610±220 








S XVI La 


4.730 


4, 


.737±0.001 


2.12±0.24 


540±320 


4.734±0.003 


1.82±0.40 


420±390 


S XVr 


5.039 


5, 


.046±0.006 


0.34±0.24 


250±170 


5.042±0.008 


0.40±0.23 


210±200 


S XVi 


5.067 


5, 


.056±0.007 


0.40±0.30 


780±440 


5.070±0.015 


0.31±0.27 


210±200 


8 XVf 


5.102 


5, 


.106±0.005 


0.44±0.15 


720±910 


5.103±0.006 


0.60±0.26 


270±260 


Si XIV L3 


5.217 


5, 


.222±0.005 


0.55±0.25 


390±510 


5.225±0.007 


0.53±0.30 


300±330 


Si XIV La 


6.183 


6, 


.192±0.001 


1.54±0.18 


550±110 


6.185±0.002 


3.43±0.47 


490±230 


Si Xlllr 


6.650 


6 


.662±0.005 


0.24±0.11 


280±310 


6.643±0.007 


0.69±0.33 


390±480 


Si Xllli 


6.669 




6.690 


<0.11 


<310 


6.705±0.010 


0.21±0.25 


160±580 


Si Xlllf 


6.740 


6, 


.746±0.005 


0.30±0.09 


800±410 


6.741±0.003 


1.12±0.34 


290±290 


Mg XII L/3 


7.106 


7 


.104±0.012 


0.36±0.13 


1450±1100 


7.109±0.007 


1.07±0.40 


690±530 


Mg XII La 


8.422 


8, 


.432±0.005 


1.40±0.33 


670±200 


8.422±0.007 


4.65±1.43 


690±370 


Mg XI(r,i,f) 


9.230 


9. 


.219±0.084 


2.31±1.01 


3940±680 
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density and excitation mecha nism (jBautista fc KallmanI 
l200Ct [Porguet fc Dubaull200C( ): each hne complex consists 
of three components, denoted r, i, f, ordered in increasing 
wavelength. At densities below the critical density we ex- 
pect R=f/i>l, and conversely. Critical densities depend 
on the nuclear charge, and range from ~ 10^ cm~^ for 
C to 10^^ cm~'^ for Fe. Gas that is photoionized is ex- 
pected to emit these lines primarily by recombination cas- 
cade, and will have G=(f-|-i)/r>l. Gas which is excited by 
electron collisions, or in which resonantly scattered emis- 
sion is important, will have G<1. The He-like complexes 
in obs. Ill all provide us with constraints on the excitation 
mechanism and gas density, although taken together they 
provide conflicting indications of these quantities. The Fe 
XXV feature peaks near 1.86 A, which is the rest energy of 
the intercombination line, but it extends beyond 1.87 A, 
the energy of the forbidden line, and there is little emission 
at 1.85 A, the energy of the resonance line. The Ca XIX 
line behaves similarly. We only measure an upper limit to 
the resonance line, but find strong intercombination and 
weak forbidden lines. However, there is excess emission at 
wavelengths greater than the laboratory wavelength of the 
forbidden line, which cannot be adequately fitted by the 
XSTAR model. Ar XVII is an intermediate case in which 
the fit accepts both a resonance line and a forbidden line. 
S XV is an intermediate case as well, though its line com- 
ponents appear less significant. Si XIII is a clearer case, 
and we have data from both observations. Figure H] shows 
the spectrum in both cases, in which the resonance and 
forbidden lines are clearly detected and the intercombina- 
tion lines are upper limits. This suggests a density below 
the critical density for the 2s'^P - 2s'^S transition, which 
is 5 xlO^"* cm^'^ for Si. The He-like line ratios from the 
heavier elements suggest that if there is a single plasma 
density it is likely near that limit. Although we cannot 
resolve Mg XI in obs. Ill, the line clearly peaks at the lo- 
cation of the intercombination line indicating an R ratio 
well below 1 and thus densities significantly above 10^'^ 
cm~^ (see also the discussion in Sect. 4.1). 



4 6 
Wavelength (Angstrom) 



Fig. 3. — Flux ratios //////// for lines measured in obs. Ill and 
II (values taken from SBII). 

The fluxes of the emission lines are shown in Table 3. 
The lines in obs. IV appear bright with fluxes between 2 
xlO^^ and 3.4 xlO~^ ph cm^^ s^^, in obs. HI ranging 
from 1.5 xlO~^ ph cm~^ s~^ to 2.1 xlO~^ ph cm"^ s^^. 

^http:/ /heasarc. gsfc.nasa.gov/docs/software/xstar/xstar. html 



The line fluxes between the two observations are not sub- 
stantially different from each other, most pre-zero phase 
fluxes seem slightly lower. This is different from what we 
observed during the high flux states (SBII) where the line 
fluxes differed substantially. The flux difference between 
pre- and post-zero phase then also was an order of mag- 
nitude larger than what we observe now. Another issue 
arises when we compare the pre-zero phase lines between 
obs. II and III (see Figure [3]). The reason we use obser- 
vation II is because here we expect the least perturbing 
effects of pileup on the lines. The source flux in obs. II 
was determined to be 4.8 x 10~^ erg cm~^ s~^ which leads 
to a continuum flux ratio of about 150. In case that we 
observe a similar photoionization regime, we would expect 
that the line ratios are very similar to this flux ratio. The 
ratios in Figure [3] appear much smaller than this expecta- 
tion. Some of this effect is certainly caused by pileup in 
obs. II which likely causes the slope from large to small 
wavelengths. Even with some of the line fluxes reduced 
by pileup effects, the line ratios fall still well short of the 
expectation from the continuum ratios, which besides the 
effects of the now reduced source luminosity we attribute 
to additional absorber activity as well as a more dramatic 
change in ionization balance. 




Fig. 4.— The Si XIII triplet from Obs. Ill (left) and IV (right). 
In both triplets the forbidden line is the strongestxi and the inter- 
combination line is weak. The slight redshift for obs. Ill is within 
the systematic uncertainty of the wavelength scale. 

3.2. Photoionization Fits with XSTAR 

The spectral analysis based on Gaussian fits indicates 
complex ionization structures. First we investigate pho- 
toionization in observation HI. For the XSTAR fits we use 
the PHOTEMIS analytic model in XSPEC, which is avail- 
able from the XSTAR website. This model generates 
a line spectrum characterized by element abundances, ve- 
locity turbulence, ionization parameter and normalization. 
Density is not a free parameter in the fits, but it can be var- 
ied by importing XSTAR model runs calculated at various 
density values. Temperature is also not a free parameter in 
these fits, since it is calculated self-consistently by XSTAR 
under the assumption that heating due to photoionization 
is balanced by the radiative cooling of the gas. The deter- 
mination of the ionization parameter is of the most interest 
in this analysis. 

The model spectrum produces all the lines listed in Ta- 
ble 3 and reproduced the listed fluxes within the quoted 
errors for over 75% of the lines. The lines are slightly 
broadened (see above) which was reflected in a turbulence 
velocity parameter of 300 km s~^ . The temperature of the 
photoionized plasma is largely set by the detection of ra- 
diative recombination continua. Since there are no signif- 
icant detections of these, the temperature of the plasma 
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TABLE 4: THE FIT RESULTS USING THE PHOTOIONIZATION (PH) AND WARMABSORBER (WA) 

MODEL IN XSTAR 

Obsid Model log N™'" b^N™^^™ I T i ^ ^ 

[cm— 2] [cm— 2] [ergs cm s^-"^] [10^® erg s^-'^ ] 



6148 PH '22.9Zoo ~ 3.0 l-58Iooo 0-001 1-219 0.79 
5478 WA (quics.) 22.3lo-} 23.9l"o ^-^-ot 0.38t°-ig -0.0077l°o°47 1.51 1.4 
5478 WA (flar.) 23.4tg:^ 22.8j:g:} '^■'^-o.l 2.62±g:}^ -0-0019±q:oo21 8-31 9^ 



a) xlO"^" erg cm"^ between 2-10 keV 

b) for a distance of 6 kpc 



must be greater than 10^ K which blends these features 
into the continuum; the temperatures of the XSTAR mod- 
els at the best-fit ionization parameter are consistent with 
this constraint. The density of the plasma is determined by 
the He-like triplets. As indicated in Sect. 13. ll the triplets 
provide conflicting trends. The fit allows a lower den- 
sity solution as suggested predominately by the Si XIII 
triplet with ^ 10^^ cm~'^, which is only in conflict with 
the density constraints from Ca XIX and Fe XXV, which 
seem to suggest > 10^^ cm~'^. An ionization parameter 
— Lx/{nr^ j) of log^=3.0±0.1 gives the best result and 
provides a cold column of 7.9±0.1 x 10^^ cm~^, a power 
law index of F = 1.58 ± 0.1, and a total X-ray flux of 
1.22x10"^° erg cm~^ . The fitted flux and power law 
index are consistent with our continuum analysis. XSTAR 
also allows some constraint on abundance and although 
some slightly higher abundances are seen for Si and Ca, 
the fit clearly required significantly enhanced Fe. 

3.3. Line Emissivities and Predictions 

It is straightforward to explore the implications of our 
measurements of emission line fluxes for models for the 
accretion flow in Cir X-1. For example, from the Gaus- 
sian fit to the lines we find a net flux of approximately 7.9 
xlO"'^^ erg cm^^ s^^ in the lines making up the Fe XXV 
K complex. At a distance of 6 kpc this corresponds to a 
line luminosity of 3.2 xlO'^^ erg s"^. If this were emitted 
by recombination in a region with uniform conditions, the 
luminosity would be: 



n^Va 



2p£lineXFeXXViyfe 



(1) 



where n'^V is the emission measure, a2p is the effective re- 
combination rate into the 2p levels of Fe XXV, eii„e is the 
line energy xp^xxvi is the ion fraction of Fe XXVI and yf^ 
is the iron elemental abundance. From XSTAR models we 
have a2pXfeXXVi = 1-95 x 10~^^ cm'^ s~^. Combined with 
the measured line luminosity, we find that the emission 
measure at log(^)=3 must therefore be ri^V = 4.9 x 10^* 
cm~^. 

This quantity can be compared with what is expected 
from an X-ray heated accretion disk coronae (ADC) in low 
mass X-ray binaries. Simple estimates for the emission 
measure distribution from an ADC are complicated, ow- 
ing to the fact that the coronal temperature, density and 
ionization are coupled by the constraints of hydrostatic 
equilibrium, thermal balance, and the radiative transfer 
of photons into the corona. Calculations for the ADC 
structure including these effects have been published by 



[London fl982^ :lBegelman et all (119831) :lKallman &: Whitd 
(198 9): ,Ostriker et all (I1991D: iKo &: KallmanI ([19941) : 



R ozaiiska et al. ( 2002D : Dimenez-Garate et al. ( 2003) . aiid 
many others. A simple estimate for the distribution of 
emission measure with ionization parameter 



dEM 



n^2TrRdR^ 



(2) 



can be obtained using an isothermal approximation to the 
disk vertical structure: n(z) = noexp{—{z/zs)'^), where 
the scale height is Zs = y^{kTR^/{GMmH)). Then: 



dEM 



(- 



kT 



1/2 



= 2.5xl0°^C^-'i38-Ri7 



65^-3^2 p-3/2jnl/2 



di ^ \GMmHRl 

where Las is the source luminosity in units of 10'^* erg/s, 
Rif is the disk inner radius in units of 10^ cm and T/^ is 
the corona temperature in units of 10^ K. 




12 3 4 

log(ionization parameter) 

Fig. 5. — Calculation of emission measure in dependence of ion- 
ization parameter for a static corona for two different luminosities, 
L = 10''^ erg (red) and L = 10^^ erg (green). 

This result is quite accurate when compared with the re- 
sults of numerical solutions for a single-stream treatment 
of the transfer of the X-rays, and for a disk around a 1 
AIq compact object extending from 10^ to 10^"'^ cm. Fig- 
ure [5] shows such a calculation, for ^33=1 (red curve) and 
7^38=0.1 (green curve). This shows that at log(^)=3, an 
emission measure that we infer from the observations can 
only be produced by an ADC if the source luminosity is 
well above 10'^^ erg/s. Although this estimate of emission 
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measure is idealized, and assumes a static corona above a 
thin disk, it is clear that the continuum luminosity needed 
to produce the line emission we observed exceeds the lu- 
minosity we infer from continuum fits by a factor almost 
two orders of magnitude. One possible explanation for this 
is that the true luminosity intrinsic to the source is much 
greater than what we observe, and that the apparent low 
state is a consequence of obscuration or collimation of the 
radiation into a direction away from our line of sight. 

3.4. The warm absorber in observation IV 

Figure [H] clearly indicates that the hard spectrum in 
obs. IV has not only changed in its continuum properties, 
but also exhibits complex absorber activity. Encouraged 
by the similarity of the line spectra above about 4 A in 
obs. Ill and IV we use the results from Sect. 13.21 in the 
modeling of the post-zero phase spectra. We adopt two 
XSTAR components to model the spectrum in obs. IV, 
the warm absorber model (WARMABS) plus the pho- 
toemission model (PHOTEMIS) determined from obser- 
vation III. The data quality is not sufficient to constrain 
all the parameters of both components unambiguously, so 
in order to reduce the number of free parameters we fix 
the photoemission component to the one fitted in observa- 
tion III. This choice also makes sense under the presump- 
tion that dips are generated by blocking the central X-ray 
source and thus exposing the accretion disk coronal spec- 
trum as o bserved in, for example, dipp ing sources like EXO 
0748-676 (jJimenez-Garate et al.|[2Q03[ ). The foremost pho- 
toionization parameters we freeze from the above analysis 
arc the ionization parameter for the photoemission, the 
normalization and the turbulence velocity parameter. In 
a second step we recognize that the obs. IV shows signif- 
icantly more variability than obs. III. We then separated 
the quiescent segments in the light curve of obs. IV (see 
Figure [T]) and the flaring segments. 




Wavelength (A) 

Fig. 6. — The total spectrum below about 5 A from obs. IV show- 
ing photoemission of the soft side but complex photoabsorption on 
the hard side. 

In the following we denominate the two levels accord- 
ingly as "quiescent" and "flaring" states. Table 4 shows 
the results of the XSTAR flts to the quiescent and flaring 
state. The model fits well to both states and produces sig- 
nificantly different results. During the quiescent state the 
warm absorber model results in a low ionization parameter 



of log ^ = 1.6 ± 0.3. The absorption feature at Fe K then 
corresponds to absorption from Fe XVI - XVIII. The cold 
column seems reduced with respect to the one observed in 
obs. Ill, but a warm column of 7.9x10^^ cm~^ appears. 
In the fiaring state the cold column is reduced as well with 
a lower warm column of 6.3 x 10^^ cm~^ but a high ioniza- 
tion parameter of log ^ = 2.7 ± 0.2. The absorber near Fe 
K has its strongest contributions now from Fe XXI-XXIII. 
In order to obtain the best fit a blueshift of approximately 
600 km s~^ is adopted. However, this value yields a flt 
which is only marginally improved (A^^ —3) compared 
with no blueshift. Similarly, tests of the fits to the obs. 
IV data with and without the emission component show 
marginal improvement (A^^ =8 for the flaring state and 
Ax^ —& for the quiescent state). Table 4 allows lists the 
total luminosities for the two WA cases. 

Figure [7] shows the dependence of the depths of ob- 
served absorption lines on ionization parameter for the 
flaring state. The fiducial optical depth is chosen arbi- 
trarily, i.e. the column density of the absorber. Plotted is 
the log of the ratio of the optical depth produced by an XS- 
TAR model wit h total column dens ity lO'^^ cm^'^ and so- 
lar abundances (|Grevesse fc Sauvai| [l999) to the observed 
depths. Most curves intersect near log ^ ^ 2.7. 

These were determined using a simple Gaussian fit to 
the lines, and the plotted error bars refiect the statistical 
error associated with these fits. This shows the range of ^ 
which is approximately consistent with the observations, 
and also illustrates the nature of the constraints: the Fe 
XXVI line constrains ionization parameter on the low side, 
and the Si and S lines constrain the ionization parameter 
on the high side. 

S I ■ ' ■ ■ I helXXVI 

Fe_XXV. 

_ Fe_XIX 



Ca_XX 

Ca_XIX 




log(l;) 

Fig. 7. — The dependence of the depths of observed absorption 
lines on ionization parameter for the XSTAR fits in the flaring state. 



4. DISCUSSION 

The at times ultraluminous X-ray source Circinus X- 
1 has been observed by nearly all X-ray missions over the 
last 30 years, and it seems that the source continues to 
provide a plethora of new features in its X-ray emission. 
The light curves, as expected from the monitoring with the 
ASM, show much less activity before but also after zero 
orbital phase. This confirms that the source has moved 
into a novel physical state. In fact, we suspect that Circi- 
nus X-1 was still in transition because since January 2007 
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until the submission of this paper the X-ray flux dropped 
under the detection hmit of the ASM. Whether the source 
is in the process of turning off or reaching a defined mini- 
mum is unknown. We sus pect the latter if we follo w the 30 
yr hght curve compiled by [Parkinson et al] (|2003f ). though 
this light curve does not provide fully continuous coverage 
during its previous low-flux periods that lasted between 
1971 and 1986. The analysis in this paper aims to identify 
the physical state of the X-ray emitting plasma and draw 
conclusions about the source's condition while it transits 
towards lower fluxes. We also put these flndings into con- 
text with the sources' previous high-flux behavior. 

4.1. Some General Remarks 

There are still many outstanding questions concern- 
ing the detailed nature of the system. B esides its iden- 
tifica tion as an accreting neutron star (jTennant et al.l 
Il986f ). there is circums tantial evidence that Cir X-1 could 
be a Z-type LMXB (ISchulz Hasinger fc Trumped Il989t 
iHasinger fc van der Klislll989l ) based on its s pectral varia- 
tion p attern in the color-color diagram (e.g.. IShirev et al.l 
119991: iDing. Qu fc Lil [20061) . However, given recent find- 
ings, it has to be a quite peculiar one for several rea- 
sons. Its long-term light curve fits more the descrip- 
tion of an X-ray transient with excursions from 1-2 
Crab to less than a mCrab (J'axkinsori et a,l, 2003) . Our 
observation s, as well as recent SWIFT and RXTE ob- 
servations ()Jonker et al.l l2007a[ ). show that Cir X-1 has 
gone from an all-time high during the 1990s into an ex- 
tended low-fiux state where the source dropped to a flux 
of 3xl0~^^ erg s~^ corresponding to a luminosity that is 
below lO'^^ erg s^^ . LMXBs that low in luminosity are 
rare and the detection of an LMXB close to the Galac- 
tic Center with a luminosit y of only 4x10*^^ erg s~^ but 
high-luminosity outbursts jMuno et al.l 12005^ may indi- 
cate the existence of a peculiar class of LMXB transients. 
Conventional Z-source do not behave that way, and even 
atoll sources are brighter in their hard states. A most 
compelling indication that Cir X-1 is indeed a neutron 
star is a recent detection of twin kHz qu asiperiodic oscilla- 
tions (QPOs) during the high-flux state (jBoutloukos et al.l 
[200^ . which seem in good agreement with Z-sources, but 
not with black-hole binaries. Furthermore, the discussion 
about the companion has r ecei ved attention again more 
recently as iJonker et al.l (|2007bl ) produced some evidence 
that the companion star of Cir X-1 is of intermediate mass 
and of types B5 to AO. This is unusual for a Z-source, as 
th e ones known all have low-mass companions. 

iJonker et al.l (|2007b[ ) constrain the companion mass to 
< 10 Mq and the angle of inclination to be > 14°. The 
angle of inclination has been under dispute s ince evidence 
for ult rarelativistic radio jets was presented (jFender et al.l 
l2004af) . An angle of < 5° became a natural choice under 
the assumption that the energy injection happened during 
the last zero-phase outburst and that there is no bend- 
ing. However, there is no evidence fo r this correlation . 
Neither recent claims of soft exces ses (llaria et al. 
nor unabsorbed soft components fPing. Qu fc Li 



2005 ) 



20061) 



In lDing. Qu fc Lil (|2006l ) it was suggested that this excess 
could come from the accretion hot spot and it large emis- 
sion volume would then rule out a large inclination. How- 
ever the spectral analysis presented seems clearly tainted 
by an unrealistically high Fe K line flux. Besides the fact 
that the continuum above the Fe K edge cannot physically 
produce such a line strength, we resolved Fe K lines in the 
high state and found them to be two orders of magnitude 
lower in strength (SBII). This leaves us with the unfor- 
tunate reality that we are still quite uncertain about the 
angle of inclination. 

If the result that the companion is of intermediate mass 
holds up, our analysis can very much rule out a possi- 
ble later B-type stellar wind as the region of origin for 
the X-ray lines. Later type B-stars with surface temper- 
atures above 10^ K can still maintain substanti al winds. 
The size of the system is about 5x10^^ cm (Ta uris et al.l 
Il999l ) which for the measured ionization parameter results 
in densities of the order of 4x10^ cm~'^, which even for 
a late B-types seems quite low. The He-triplets suggest 
much higher densities than that which places the ionizia- 
tion region much closer to the X-ray source. We also mea- 
sure signiflcantly enhanced Fe in the photoemission, which 
is quite hard to argue for in an early type stellar wind. 




can be substantiated in our high-resolution spectra dur- 
ing the low-flux sta te here and in the hig h-flux state years 
ago (see also SBII. [Callowav et al.11200 7'). and any claims 
referring to the angle of inclination likely do not apply. 



Fig. 8. — The range of R and G ratios for the Ho-hkc Une emission 
in Obs. III. 

Before we discuss emission regions in more detail, we 
have to address our difficulties determining a plasma den- 
sity from the He-like triplets. Though we are somewhat 
plagued by bad statistics and some uncertain systematic 
effects, we still expected a more consistent picture. To date 
our ability to diagnose the high-density regime has always 
by limited by the unavailability of data and these observa- 
tions are a rare exception. The Mg IX triplet indicates an 
R-ratio << 1 which puts a lower limit on the density of 
lO^'^ cm"'^. Figure E] shows a cartoon relating the R-ratios 
with the corresponding G-ratios. High-density regimes are 
a consistent solution for all triplets but Ar XVII, which is 
however statistically weak. Ca XIX, on the other hand, 
though it would be consistent with high density, does not 
fit into the recombining regime and appears more coUi- 
sional. Fe XXV provides the tightest solution for higher 
density, but is also close to a coUisional solution. Though 
all in all a density of ~ 5 x lO^'' cm~^ in a recombining 
gas seems most plausible, we want to caution that there 
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are remaining inconsistencies, which, given the turbulent 
nature of the system, could indicate time-dependent as- 
pects and a situation where our equilibrium assumption 
may not be entirely valid. 

4.2. The Photoionized Disk Corona 

Quite evidently, the P Cygni lines found during the early 
Chandra observing cycles (BSI, SBII) arc not present in 
the new spectra. So far these line profiles still stand unique 
for this class of objects. The last time we observed these 
P Cygni lines was during a short observ ation of Cir X- 

1 at phase 0.25 during the year 2001 (jGallowav et al.l 

l2007f ). The interpretation at the time was that the profiles 
are a result of a preferentially equatorial wind driven by 
the central source. In observation III, the overall source 
flux was reduced by over two orders of magnitude with re- 
spect to the 2000 observations and one of the questions 
to answer is if this is a true reduction in source lumi- 
nosity or an effect of increased circumstellar/circumdisk 
absorber activity. Our direct comparison of the ratio of 
the line fluxes with the ratio of the continuum fluxes sug- 
gests that absorption is playing an important role, even 
though that comparison is tainted by the pileup problems 
during the high-state observations. On the other hand, 
photoionized plasmas are complex and we have to look 
at the change in ionization balance between the high and 
now low flux phases represented by the ionization param- 
eters. Although we were not able to directly measure ^ in 
obs. I and II we could in fer parameters well in excess of 10"' 
using calculations from ' Kallman fc Bautistal (|200l[ ). Our 
fits to obs. Ill resulted in an ionization parameter of 10"^, 
which is an order of magnitude lower. Thus significant 
dimming of the X-ray flux thereby changing the ioniza- 
tion balance provides a solution that would not require 
significant additional absorption. Thus the comparison of 
line and continuum fluxes and the comparison of ioniza- 
tion balances have quite different implications. In that 
sense the observed X-ray emission may not provide a full 
account of the true energy balance near the central source 
and we have to acknowledge the fact that there has been 
no significant long-term increase in the observed cold col- 
umn between the high and the low flux states for the years 
2000 and 2005, respectively. 

At this point it seems necessary to invoke some model 
assumptions about accretion dynamics. In order to test 
the assumption that the photoemission region could be 
the accretion disk we compare results from our XSTAR 
fits of obs. Ill to some simple model predictions. As a ba- 
sis for these tests we use som e of the models described by 
iJimenez-Garate et al.l (j2002l) . who modeled an accretion- 
disk atmosphere and corona photoionized by a central X- 
ray continuum source and also discuss feedback mecha- 
nisms between disk structure and the central illumina- 
tion. In t his model we assume an o ptically thick stan- 
dard disk (jShakura &: SunvaevI I1973D with a heated at- 
mosphere, which is separated from the disk by a r = 1 
surface, sometimes called photosphere. The atmosphere is 
strongly heated resulting in evaporation into a hot corona 
and eventually leads to outfiows. The plasma above the 
photoshere is heated up to 10^ K and cools through atomic 
line emission. The most important parameter to investi- 
gate the state of the plasma is the ionization parameter ^. 



In order to determine the source luminosity we assume 
a source distance of 6 kpc. The distance is not well known 
but it is widely agree d that the sour ce is not closer than 
this. The claim bv ilaria et al. (l2005l) of 4.1 kpc has been 
recently dismissed by iJonker et al.l ( 2007bD who actually 
favor a distance of ~8 kpc. For the density n we have 
some indication from the He-like triplets, but plausibility 
arguments are still needed. 

In the inner disk under the idealistic assumption that 
viscous heating and radiative heating from an illuminating 
source is locally expressed by a blackbody, radiative heat- 
ing should domin ate over viscous he ating at radii larger 
than 2.3x10* cm ()Vrtilek et aLlligOO V With a luminosity 
of 8x10'^^ erg s~^ we expect an upper limit to the density 
of ~ 10^^ cm~^. The strength of the forbidden line in the 
Si XIII argues against such a high density and if we employ 
our projection of ~ 5 x lO'^'* cm^'^, we place the emission 
region at a radius of ~ 10^ cm. During the high flux in the 
year 2000 we projected a launching radius > 10^*^ cm and 
densities > 5 x 10^^ cm~^, which given the higher source 
luminosity seems consistent. 

As shown in Sect. 13.31 luminosities which can heat an 
ADC have been observed during the high-state of the 
source. The luminosities inferred during the observations 
reported here are signiflcantly lower. This suggests that 
either there are problems with the static corona approx- 
imation or that there are unobserved photons or other 
heating mechanisms responsible for the coronal excitation. 
We do not see any obvious signs for winds in our data like 
in the year 2000 observations, however there still seems to 
be jet activity in the radio (P. Jonker, priv. comm.). On 
the other hand, the problem of insufficient heating power 
for ADCs as derived from X-rays does not seem to be 
unique to Cir X-1 . The emission measures determine d in 
classi cal ADC sources such as 4U 1822-37 ((Cotta m et al.l 
1200 ID . EXO 748-676 (iJimenez-Garate et al.l 12003 



20031) and HerX- 



2S 0921-63 (iKallman et al.l 

1 (|Jimenez-Garate et al.ll2005fl are greater than the ADC 
model can produce based on the observed luminos- 
ity. In the case of Her X-1 it has been suggested that 
the energy range is indeed much dom inated by absorp- 
tion (|Kuster et al.ll2005l : iJi et al.|[2007l) and the luminosity 
in its low-state may be higher than observed. If this is 
the case, then the ionization balance we observe in the 
low-state of Cir X-1 should also be largely dominated by 
obscuration rather than a true reduction in source lumi- 
nosity and as such would suggest a genuine analogy to a 
Seyfert II. 

All the sources mentioned above have high angles of 
inclination (85° < i < 90°) which given the strong simi- 
larity of the low-state spectral properties in Cir X-1 is a 
confirmation of the previous X-ray projection of a higher 
angle of inclination. However, there are two distinct dif- 
ferences between these sources and Cir X-1 . One is the 
fact that we have a history of the source of once being 
very luminous exhibiting powerful winds, the other one is 
its jet activity which even in the low-flux state seems to 
be prevailing. Generally almost two orders of magnitude 
in luminosity are not easily hidden by absorption and the 
relative constancy of Nh in Cir X-1 over the years is diffi- 
cult to reconcile with a simple obscuration paradigm. The 
fact that the luminosity changed from a level that seemed 
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sufRcient to sustain a static corona in the high-state to a 
gross deficiency in the low-state without much change in 
absorption properties between obs. II and III also leads 
to questioning whether the assumption that the X-ray lu- 
minosity is the dominating heating component and/or the 
assumption of a simple static corona is correct. 

Jet and wind activity indicates that the disk in Cir X- 
1 may be cooling more effciently than predicted from stan- 
dard models and additional energy is transferred into the 
corona or the jet base. The connection of compact ADCs 
and the base of jets has been ma de recently in a study by 
iMarkoff. Nowak fc Wili^ (|2005[ ). Though it is well be- 
yond the scope of this paper to go into details, we have to 
address the possibility that other mechanisms than X-ray 
ilUumination can heat the corona. Recent studies clearly 
emphasize the role of magnetic activity. 

Relativistic jets a re a fundamental aspect of accr etion 
on to black holes (jFender. Belloni fc Galiol l2004bD and 
much emphasis in recent years has been devoted to inves- 
tiga te connections bet w een compact ADCs and jet activ- 
ity iMerloni fc Fabian' (2002^; 'Markoff. Nowak fc Wilmi 
(|2005[ ): iKrohk. Hawlc y fc Hirose (2007.) . The for- 
mati on of jets may arise from radiatio n pres- 
sure (lLi2iden-Beir'1978'; 'Begelman fc ReedllSSl 

or mag- 
netic fields (Blandfordfc Znayck 1977) and since we deal 
with X-ray luminosities well below L^dd we may argue 
for the latter as the underlying scenario. Viscous turbu- 
lence in the accretion flow is now under stood to be due to 
the m agnetorotational instability (MRI, iBalbus fc HawlevI 
Il998f ) and it may as well also the source for a ma gnetically 
dominated corona (jMerloni fc NavakshinI I2006D . Obser- 
vational evidence for MRI effects in GROJ1650 -55 has 
recently been presented by iMiller et al.l ()2006cD . The 
low-flux observations presented in this paper fit most cri- 
teria that qualifies the low-flux state as some form of 
low/h ard state as seen in accretion black hole binaries 
(see iFender. Belloni fc Gallol l2004bt iMiller et all l2006bl 
and references therein). 

4.3. The Absorber Region 

The fact that there is a significant warm absorber in 
the post-zero phase spectra is obvious through the broad 
absorption feature around 1.9 A. In hot gases, where the 
electrons are generally stripped out of their atomic bond, 
absorbing fractions are usually generated by the formation 
of a significant number of He-like and H-like ion states. 
Such absorption has been observed in many accreting X- 
ray binaries, some show ionization fractions only for iron, 
some throughout the X-ray band down to oxygen. Circi- 
nus X-1 in this respect was pioneering as the Chandra first 
observations at the end of the high flux period showed 
blueshifted absorption lines for most detectable H- and 
He-like ions species. On top of this phenomenon we could 
show in SBII that the equivalent widths of the hot absorp- 
tion is tightly linked to the slope variation of the underly- 
ing continuum. 

The changes we see in the warm absorber are a bit more 
subtle but point in a similar direction. We now observe an 
entirely different ionization regime and the atomic physics 
is more complex due to the fact that in this absorber we do 
not have simple He-like and H-like ions. Atomic data cov- 
ering this ion regime has not been available until recently 



when lKallman et al.l (|2004D updated the XSTAR database 
specifically in the Fe K band. Model simulations indicated 
that under very specific circumstances Is-np autoionizing 
transitions become dominant. The results are Fe XVII - 
XIX absorption features near 1.9 A and the Fe K edge. 
The appearance of this edge becomes increasingly diluted 
as the ionization parameter increases. The model regime 
then predicts warm columns of 10^^'^ cm^^ and ionization 
parameters well below 100. 

The unique ionization regime that produces the ob- 
served absorption makes the interpretation of the absorp- 
tion relatively simple and straight forward. An optically 
thick cool and distant plasma is illuminated by a moder- 
ate source luminosity producing low-ionization ions. What 
makes the situation more complex is the fact that we ob- 
serve at least four absorber regimes: a cold absorber, a 
"luke-warm" absorber during quiescence, a warm absorber 
during flares, and last but not least a hot absorber during 
both phases. 

The cold absorber: In all our observations we detect 
cold photoelectric absorption at levels significantly above 
the expected line-of sight-columns. In obs. HI we observe 
about 8x10^^ cm^'^, in obs. IV it is lightly lower in qui- 
escence. Compared with obs. I and II from the high flux 
state of the year 2000 these numbers are not so much dif- 
ferent. If we compare our results from the partial coverage 
fits with the ones in SBII, but als o with the ones fro m the 
ASCA spectra a few years earlier (jBrandt et al.lll996l ) then 
cold absorption remains at a rather steadly level for long 
time periods. The cold absorber should be devoid of sig- 
nificant ionization fractions above Fe X in the X-ray band 
which limits the ioniziation parameter to below 1 and puts 
the material at a distance consistent with the outer accre- 
tion disk, i.e beyond 10^^ cm. The amount of cold material 
so far seems to be fairly invariable. Not only is it rather 
constant at zero orbital phases, but the edge optical depth 
remained constant th roughout other binary orbit phases 
(jCallowav et al.lF2007f ) indicating similar amounts of cold 
material at all times. 

The luke-warm absorber: The fits at quiescence 
shows an enormous warm column of 8x10^^ cm^^. This 
warm absorber shows up in addition to the cold column. 
The slight reduction in cold column can be easily seen at 
the expense of a slightly increased source luminosity. It is 
evident, however, that this large mount of material has not 
been in the line of sight during the pre-zero phase obs. HI 
as otherwise it could have only been exposed through a 
luminosity decrease, which would have rendered the warm 
absorber to a cold one. There are two possible time lines 
of how the matter got there. Either the material crept 
up over the time period between the observations, or we 
see newly accreted matter just after the incoming neutron 
star-disk system re-connected with the companion. Given 
the invariability of the cold column over several observ- 
ing periods and that a sudden inc rease in mass accretion 
at zero orbital phase is expected (jTauris et al.l [1999) we 
consider the latter possibility as more likely. During the 
high fiux phase we did not see this effect as there the ma- 
terial was almost completely ionized at all times. This 
material, on the other hand, is then likely related to the 
materi al causing heavy dipping activity during that flux 
period iBrandt et all (|1996f ): IShirev et all (119990 . 
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The warm absorber: The reason why we separate 
the absorber in quiescence and during flaring periods is 
that due to the large increase in source luminosity the 
absorber also experiences a change in ionization balance 
which also create higher transition Fe ions. The amount 
of warm column is now only 6.3x10^^ cm~^ which is an 
order of magnitude lower than during quiescence. There 
may two causes we can consider for the reduction. First, 
the source is actually accreting and thus mass moves to- 
wards the central source. This flow is usually subsonic 
(|Frank. King fc RaindfTOQ^) and matter moves out of the 
ionization regime on timescales much larger than the entire 
observing time. The light curves in Figure[T]feature flaring 
periods of only a few hours. More likely is a feedback effect 
from direct accretion. Accreted matter is converted into 
X-rays further heating the absorber. The increase in cold 
absorption during flaring may simply reflect the change in 
ionization state of accreted matter. 

The hot absorber: The fact that we also detect H- 
and He-like absorbers in the spectra also indicates the ex- 
istence of a hot absorption component which reveals itself 
in resonance absorption of Ca XX, Fe XXV, and XXVI. 
The positions of the lines are not blue-shifted as observed 
during the high flux phases indicating that this absorber is 
not a wind but a stationary component. The hot absorber 
is likely part of the warm absorber reflecting the upper 
end of the spectrum of ionization fractions. 

The light curve of obs. IV shows that final accretion 
has two components. In quiescence matter is supplied at 
a fairly steady level. During flares accretion also comes 
in chunks of different sizes which is reflected by the very 
short flare onsets and types. Average rise times are less 
than a few minutes. Large flares are accompanied by a 
staccato of many subflares indicating rather turbulent ac- 



cretion. Most remarkable is the correlated change of ab- 
sorber column and ionization fraction with the luminosity. 
Although our data only allow us to analyze the integrated 
spectra of quiescent and flaring intervals, this analysis in- 
dicates that all the absorbers are connected and the source 
behaviour may be explained by varying illumination levels 
of the absorbers combined with varying obscuration levels 
of the central source. 

There are many places in the system where this can hap- 
pen. Based on the observed ionization fractions, most of 
the material should be at larger radii from the source than 
the photo-emission regions. This is likely also true for the 
hot absorber, which can be seen by the fact that during 
flares, the flux of most emission lines are reduced. The 
increase in source luminosity cannot remove these ioniza- 
tion fractions which should survive at even much higher 
luminosities as we have observed during obs. II (SBII). 
Although we cannot rule out other more exotic options a 
most logical place for the absorbers is the mid-and outer 
disk atmosphere/corona which in turbulent systems such 
as Cir X-1 can puff up to fairly large scale heights. Even 
in this case, though, we would limit incl ination to higher 
than 60° (see also iJimenez-Garate et a l. (2002)) as we 
proposed perviouslv (BSI. SBII. ISchulz et al.n2006<) . The 
variable absorber activity we observe together with the 
complex play of warm column and luminosity (see Table 4) 
again underlines the importance of obscuration effects and 
the Seyfert II analogy. 

We thank all the members of the Chandra team for 
their enormous efforts. We gratefully acknowledge the fi- 
nancial support of CXC grant GO0-1041X (WNB, NSS) 
and Smithsonian Astrophysical Observatory contract SVl- 
61010 for the CXC (NSS). 
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